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ABSTRACT 

There exists a significant population of broad line, z ~ 2 QSOs which have heavily 
absorbed X-ray spectra. Follow up observations in the submillimetre show that 
these QSOs are embedded in ultraluminous starburst galaxies, unlike most unab- 
sorbed QSOs at the same redshifts and luminosities. Here we present X-ray spectra 
from XMM-Newton for a sample of 5 such X-ray absorbed QSOs that have been 
detected at submillimetre wavelengths. We also present spectra in the restframe 
ultraviolet from ground based telescopes. All 5 QSOs are found to exhibit strong 
C IV absorption lines in their ultraviolet spectra with equivalent width > 5A. The 
X-ray spectra are inconsistent with the hypothesis that these objects show nor- 
mal QSO continua absorbed by low-ionization gas. Instead, the spectra can be 
modelled successfully with ionized absorbers, or with cold absorbers if they posess 
unusually flat X-ray continuum shapes and unusual optical to X-ray spectral en- 
ergy distributions. We show that the ionized absorber model provides the simplest, 
most self-consistent explanation for their observed properties. We estimate that 
the fraction of radiated power that is converted into kinetic luminosity of the out- 
flowing winds is typically ~ 4 per cent, in agreement with recent estimates for the 
kinetic feedback from QSOs required to produce the M — a relation, and consis- 
tent with the hypothesis that the X-ray absorbed QSOs represent the transition 
phase between obscured accretion and the luminous QSO phase in the evolution 
of massive galaxies. 
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1 INTRODUCTION 

The prevalence of black holes in present-day galaxy bulges, 
and t he proportionality betw e en black hole and spheroid 
mass (|Magorrian et al.. 19981 : iMerritt fc Ferrarese. 2 001') 
implies that the formation of the stars in a present- 
day galaxy spheroid was fundamentally connected to the 
growth of the black hole as an active galactic nucleus 
(AGN). In fact, an evolutionary connection between the 
star-forming and AGN -dominated phases o f galaxies has 
long been posited (e.g. ISanders et al.. 19881 ). 

It is now known that a significant fraction of cosmic 
star formation took place in submillimetre galaxies 
(SMGs): dusty, ultralumi nous, and highly obscured 
galaxies at high r edshift ( Sm all. Ivison fc Blain. 19971 : 
iHughes et al.. 1998h . Although some distant QSOs 
have been detected as sub millimetre sources in 
pointed observa tions (e.g. iMcMahon et al.~ 1999: 
IPage et al.. 200 J l. the overlap between the QSO and 



SMG populations is re latively small jFabian et al.. 200C 
ISevergnini et al.. 2000 



Almaini et al.. 2003 



Iwaskett et al.. 20031 : iLutz et al.. 20ld ). In an attempt 
to uncover the relat ionship between q uasi-stellar objects 
(QSOs) and SMGs, IPage et al. f2004h observed matched 
samples of X-ray absorbed and unabsorbed QSOs at 
850/im with SCUBA. These observations revealed a 
remarkable dichotomy in the submillimetre properties 
of these two groups of sources: a substantial fraction 
(~ 50 per cent) of X-ray absorbed QSOs are submil- 
limetre galaxies, while X-ray unabsorbed (i.e. normal) 
QSOs are not. This suggests that the two types are 
linked by an evolutionary sequence, in which the X-ray 
absorbed QSOs correspond to the transition between 
the main star-for ming phase and t he QSO phase of a 
massive galaxy (|Page et al.. 20041 : IStevens et al.. 20051 : 
ICoppin et al.. 20081 ). The space density and luminosities 
of the X-ray absorbed QSOs indicate that this transitional 
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phase is relatively brief, ~15 per cent of the duration of 
the luminous QSO phase. 

However, the physical driver for this transition be- 
tween obscured star-forming galaxy and luminous QSO 
is unknown. Furthermore, the nature of the X-ray ab- 
sorption in the X-ray absorbed QSOs remains puzzling. 
These objects are characterised by hard, absorbed X-ray 
spectra, but they have optical/UV spectra which are typ- 
ical for QSOs, with broad emission lines and blue con- 
tinua. Assuming that their hard X-ray spectral shapes 
result from photoelectric absorption from cold material 
with solar abundances, the column densities are of order 
10 22 cm -2 . The X-ray absorption could be due to gas lo- 
cated within the AGN structure, or from more distant 
material in the host galaxy. Wherever the absorber is lo- 
cated, it appears to contain very little dust: for a Galactic 
gas/dust ratio, the restframe ultraviolet spectra would be 
heavily attenuated by such large columns of material. 

Therefore in order to investigate the X-ray ab- 
sorption in these objects, we have obtained deep (50- 
lOOks) XMM-Newton observations of five X-ray ab- 
sorbed Q SOs from our sample of hard -spectrum Rosat 
sources (|Page. Mittaz fc Carrera. 200lh that were de- 
tected as powerful submillimetre sou rces with SCUBA 
(|Page et al.. 200 ll ; IStevens et al.. 2005h . 

Throughout this paper we assume Ho = 



70 km s" 1 Mpc 



0.3, and Qa 



0.7. We 



define a power law spectrum such that /„ tx v~ a \ we 
assign the energy index a the subscript X when it refers 
to the X-ray spectral slope, O when it refers to the 
optical/ultraviolet spectral slope and OX when it refers 
to the optical to X-ray spectral slope. 
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Figure 1. The optical spectra of RXJ005734, RXJ094144, 
RXJ121803, RXJ124913 and RXJ163308. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical spectra 

RXJ005734 was observed on the 22nd November 1998, 
with the ESO 3.6m telescope at La Silla, Chile. The 
EFOSC 2 spectrograph was used with the 5000A blaze, 
300 lines mm -1 grism covering the wavelength range 
3800A to 8000A, with a spectral resolution of 20A FWHM. 

RXJ094144, RXJ121803 and RXJ163308 were ob- 
served on the 3rd-5th March 1998 with the William Her- 
schel Telescope at the Observatorio del Roque de los 
Muchachos, La Palma. The dual arm ISIS spectrograph 
was used with the 5400 A dichroic, the R158R grating on 
the red arm and the R300B grating on the blue arm. The 
two ISIS arms combined give continuous coverage of the 
wavelength range 3200 A to 8100 A with resolution of 5 A 
FWHM in the blue and 9 A FWHM in the red. 

In both observing runs, arc-lamp spectra and observa- 
tions of spectrophotometric standard stars were taken for 
wavelength and flux calibration. The spectra were reduced 
and calibrated using standard routines in IRAF. 

RXJ124913 was not observed because it is a 
well known broad absorption line QSO (BALQSO) 
with optical spectra available in the literature; 
the spectrum shown in Fig. [1] is taken from 
lJunkkarinen. Burbidge fc Smith ("19871 ). 



2.2 X-ray spectra and ultraviolet images 

The XMM-Newton observations used in this analysis are 
listed in Table Q] The data were reduced with XMM- 
Newton SAS version 6.5, and the calibration release of 
June 2006. After initial processing of the EPIC MOS 
and pn data, spectra were extracted from a circular re- 
gion of 15 arcsecond radius, centred on the target. Back- 
ground was obtained from a circular region of radius 
150 arcseconds, within which all significant sources were 
masked out. Appropriate response matrices and effec- 
tive area files were generated with the SAS tasks rmf- 
GEN and ARFGEN. For each object, spectra and response 
matrices were combined u sing the method described in 
IPaee. Davis fc Salvi (20031 ). To determine how the fluxes 
of the sources compared to those obtained in our original 
Rosat survey, we fitted the spectra below 2 keV with a 
powerlaw and fixed Galactic absorption. For RXJ005734, 
RXJ121803 and RXJ163308 we obtain 0.5 - 2 keV fluxes 
of 1.2 x 10~ 14 , 1.5 x 10" 14 and 1.8 x 10~ 14 erg cm" 2 s" 1 
respectively, consistent with t hose obtained with Rosat 
l|Page. Mittaz fc Carrera. 200Cj ). For RXJ094144 we ob- 
tain a 0.5 — 2 keV flux of 5 x 10" 15 erg cm -2 s" 1 , ap- 
proximately a quarter of its flux during our Rosat sur- 
vey, while for RXJ124913 we obtain a 0.5 - 2 keV flux of 

1.3 x 10~ 14 erg cm" 2 s" 1 , approximately half the Rosat 
flux reported bv lPage. Mittaz fc Carrera" (20001 ). 

Simultaneous with the X-ray observations, the XMM- 
Newton Optical Monitor (XMM-OM) took deep images 
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Table 1. XMM-Newton observations. Exposure times refer to the clean exposure times after periods of high background have been 
excluded. The UVW1 magnitudes for RXJ094144 and RXJ121803 come from the XMM-OM observations performed simultaneously 
with the X-ray observations, and are in the XMM-OM Vega system. 





Obs ID 


date 




F 11 LCI S 




Exposure time (ks) 


T IVW1 
U V VVl 








pn 




TV/TOCO 


pn 


MOS1 


MOS2 


mag 


RX J005734.78-272827.4 


0302310301 


18 Dec 2005 


thin 


thin 


thin 


40.2 


47.7 


47.3 




RX J094144.51+385434.8 


0203270101 


17 May 2004 


medium 


medium 


medium 


32.6 


40.5 


40.8 


21.63 ±0.12 


RX J121803.82+470854.6 


0203270201 


1 June 2004 


thin 


thin 


thin 


30.0 


36.9 


37.2 


20.63 ± 0.06 


RX J124913. 86-055906.2 


0060370201 


11 July 2001 


thin 


thin 


thin 


25.0 


39.6 


39.6 






0203270301 


15 July 2004 


medium 


medium 


medium 


34.3 


42.5 


42.9 




RX J163308. 57+570258.7 


0302310101 


7 Sept 2005 


thin 


thin 


thin 


5.5 


9.2 


8.8 






0302310501 


23 Oct 2005 


thin 


thin 


thin 


10.5 


19.7 


20.0 





111 the UVW1 filteifl which has an effective wavelength 
of 2910A. These observations were obtained primarily to 
provide photometric redshift constraints for objects sur- 
rounding the QSOs, and for the three highest redshift 
objects the flux in this band is severely affected by the 
Lyman break. For RXJ094144 and RXJ121803 however, 
the Lyman break cuts in only at the extreme blue of the 
UVW1 transmission, and will have a minor impact on the 
photometry. The data were reduced using the XMM-OM 
SAS version 6.5. The UVW1 magnitudes of RXJ094144 
and RXJ121803 (in the XM-OM Vega system) are given 
in Tabled] 



3 RESULTS 

3.1 Absorption in the rest- frame ultraviolet 

Broadband optical spectra of the five X-ray absorbed 
QSOs are shown in Fig. [I] It is notable that in all five 
objects, at least one absorption line is superimposed on 
the broad CIV emission line or on the continuum to the 
blue of this line. To examine this more closely, we show the 
regions around the C IV emission line in more detail in Fig. 
[2] In RXJ005734, there are two absorption lines apparent, 
with outflow velocities of 4800 and 8100 km s , assum- 
ing that they are CIV. In the spectrum of RXJ094144, 
at least one absorption line is present, apparently inflow- 
ing with a velocity of 700 km s" 1 assuming that it is 
C IV. As this spectrum is relatively noisy, we are unable 
to confirm whether the tentative absorption features to 
the blue of CIV are real or statistical fluctuations. There 
are two significant absorption features in RXJ121803, 
one inflowing at 1200 km s" 1 , and one outflowing at 
1600 km s _1 . Mgll absorption is detected from the in- 
flowing component. For RXJ 124913 we adopt an emis - 
sion line redshift of z = 2.236 (jHewitt fc Burbidee. 1991 ) . 
although the reported emis sion line redshift of this 
source ranges from % = 2.22 (|Boksenberg et al.. 197ct ) to 



1 UVW1 images were not obtained during the July 2001 ob- 
servation of RXJ124913, in which the UV grism of XMM-OM 
was used instead. 



z = 2.244 (jHill. Thompson fc Elston. 19931 ) depending on 
which emission lines the redshift is based on. Relative to 
this emission line redshift, the CIV broad absorption line 



with a FWHM 
(jBoksenberg et al.. 19781 ). In RXJ163308 



has a central ve locity of - 16700 km s 
of 5700 km s" 1 
there are four absorption lines visible, but two of these are 
likely to be due to an intervening absorption line system 
with z = 1.482, for which a number of other transitions 
are detected in the optical spectrum. The remaining two 
lines have outflow velocities of 800 and 7700 km s _1 . Both 
are accompanied by Lya absorption with consistent ve- 
locities, and absorption from Si IV is detected in the 7700 
km s _1 component. The CIV absorption line properties 
of the five QSOs are listed in Table [2] 

The five QSOs have strong rest-frame ultraviolet con- 
tinuum emission. Ultraviolet spectral slopes were mea- 
sured from the optical spectra, longward of Lya, after 
masking emission and absorption lines. The ultraviolet 
spectral slopes are listed in Table [2] With the exception of 
RXJ094144, the ultraviolet continua of the five QSOs have 
spectral slopes in the range < olq < 1, which is typical 
for optically-selected QSOs l|Francis et al.. 199ll ). indicat- 
ing little or no dust reddening intrinsic to these objects. 
On the other hand, as can be seen in Fig.[T]the rest-frame 
ultraviolet continuum of RXJ094144 is much redder than 
those of the other objects. The continuum slope in this ob- 
ject is a — 2.7, outlying from th e distribution of opt ically- 
selected QSOs at the 3<r level dFrancis et al.. 199ll "l. and 
suggesting that this object may be somewhat dust red- 
dened. The ultraviolet spectral shape can be matched by 
redde ning the median Sloan Digi tal Sky Survey QSO spec- 
trum (|Vanden Berk et al.. 200lh by E(B-V)=0.20 using 
the SMC reddening law (|Pei. 19921 ). a similar level of ex- 
tinction as observed in the nearb y starburst/BALQSO 
IRAS 07598 + 6508 (|Lipari. 19941 ). and suggesting ~1.4 
magnitudes of extinction at 2500 A in the restframe. The 
UVW1-B colour of RXJ094144 (where t he B magnitude 
comes from lPage. Mittaz fc Carrera 200ll ) is in reasonable 
agreement with this reddened template: after correction 
for Galactic extinction the UVW1-B colour of RXJ094144 
is 0.47 ± 0.18 compared to UVW1-B=0.71 predicted by 
the reddened template. In contrast, RXJ121803 has a 
much bluer UVW1— B=— 0.45 ± 0.11 after correction for 
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Table 2. Characteristics of C IV ultraviolet absorption line sys- 
tems and ultraviolet continuum slopes a a . Negative velocities 
correspond to blueshifted lines and positive velocities corre- 
spond to redshifted lines. 



QSO 


Velocity 
(km s _1 ) 


C IV FWHM 

(km s _1 ) 


CIV EW 

(A) 


a 


RXJ005734 


-4800 


500 


5 


0.6 


RXJ005734 


-8100 


500 


2 


0.6 


RXJ094144 


700 


800 


12 


2.7 


RXJ121803 


1200 


900 


15 


1.0 


RXJ121803 


-1600 


900 


14 


1.0 


RXJ124913 


-16700 


5700 


70 


1.0 


RXJ163308 


-800 


700 


11 


0.1 


RXJ163308 


-7700 


700 


10 


0.1 



Galactic extinction. For comparison, the median SDSS 
QSO template would have U VWl-B=-0.68 at the red- 
shift of RXJ121803. 



3.2 Modelling the X-ray spectra 

Analysis of the X-ray spectra was performed using SPEX 
version 2.OO.110- As a starting point, the XMM-Newton 
spectra were fitted with a similar mo del to that consid- 
ered bv lPage. Mittaz fc Carrera (2003 ): a power law with 
a fixed ax of 0.98, fixed Galactic absorption, and cold 
photoelectric absorption at the QSO redshift with Nh 
as a free parameter. The value of ax = 0.98 was cho- 
sen because it is the mean spectral index of QSOs in the 
XMM-Newton energy range and at s imilar X-ray flux lev - 
els to our targets, as determined by Mateos et al. (2 005) 
and confirmed in more recent work (jMateos et al.. 2010h . 
The results are given in Table [3] Although some of the 
individual objects are successfully fit by this model, oth- 
ers are not (RXJ121803 and RXJ124913), and the total 
X 2 /v = 225/161 for the sample of 5 objects is not accept- 
able. Therefore, we attempted to fit the spectra with two 
models, one in which the continuum spectral slope was al- 
lowed to vary, and one in which the ionization parameter 
of the absorber was allowed to vary. 

As seen in Table [3J allowing the continuum spectral 
indices or the ionization parameters of the absorbers to 
vary produced acceptable values of x* ' l v f° r the 5 objects 
individually, and as a sample. Since the two models are 
equivalent in terms of x 2 j v i we can only discriminate be- 
tween the two models by looking to see which provides the 
more physically plausible and self-consistent explanation 
for the spectra. In particular, we can test the hypothesis 
that the underlying spectra of these objects are reason- 
ably typical for QSOs, and their unusual observed spectral 
shapes are due to absorption. Two properties in particular 
are important in this context: the spectral index of the un- 
derlying X-ray spectrum, and the ratio of optical to X-ray 
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Figure 2. Close up of the region around the CIV 1550 A emis- 
sion line in our five X-ray absorbed QSOs. For comparison, 
the grey line s hows the median Sloan Di gital Sky Survey QSO 
template from IVandcn Ber k et al. (2001) 



flux, which we parameterise as the equivalent aox, which 
is the energy index of the power law that would connect 
the restframe 2500 A and 2 keV flux densities. 

Beginning with the X-ray spectral slope, the most ap- 
propriate compa rison sample of norm al QSOs is the sam- 
ple examined bv lMateos et al.' (20051 ). which is an X-ray 
selected sample with a similar flux limit to the ROSAT 
hard spectrum sample from which our targets were drawn, 
and measured in the same energy range, and with the 
same instrumentation [XMM-Newton EPIC) as our tar- 
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Table 3. Model fits to the XMM-Newton spectra. The models are a combination of a power law (PL), and either a cold photoelectric 
absorber (CA) or warm, ionized absorber (WA) at the redshift of the AGN. In addition, a zero-redshift cold absorber with column 
density fixed at the Galactic value is included in each model. Parameters marked 'f' are fixed in the fit, and '*' is used to indicate 
an uncertainty in a model parameter that is not constrained by the fit. Njj has units of cm" 2 and £ has units of erg cms" 1 . 

RXJ005734 RXJ094144 RXJ121803 RXJ124913 RXJ163308 Total 



PLxCA 
a x 

logoff 

&OX 
X 2 /» 



0.98-f 
22.0 
1.70 

41/36 



+0.1 

o.i 

+0.05 
0.06 



0.98-f 

0.2 
0.2 



22.1 



1.48 



+0.30 



0.12 

15/12 



0.98-f 
21.8 
1.32 

64/32 



+0.1 
0.1 
+0.06 
0.05 



0.98J f 
22.2 
1.89 

82/53 



+0.1 
0.1 
+0.09 

0.09 



0.98^ 
21 5+ - 2 

1 44+0-05 
^ -0.04 

22/28 



225/161 



PLxCA 

a x 

\ogN H 
aox 



n S9 +0.10 
U-»2_ og 



21.8+0-j 
, --,+0.05 
*-■' ^-0.06 



39/35 



0.57 
21.5 
1.54 
11/11 



+0.19 
-0.17 
+0.4 
0.9 
+0.52 
0.13 



+0.07 
0.06 
+21.9 

+0.07 
0.06 

29/31 



0.41 
0.00 
1.39 



0.53 
21.5 
1.96 
50/52 



+0.07 
-0.07 
+0.2 
0.3 
+0.10 
0.10 



0.94 



+0.15 



21.4 



+0.4 



1.43 



+0.05 



0.04 
22/27 



151/156 



PLxWA 
ax 
log? 
\ogN H 

aox 
X 2 /v 



0.98-f 



2.1 



+0.2 
0.3 
+ 0.2 



22 7~* 
zz -'-0.2 

l-68_ 06 
34/35 



0.98* 
7+0.4 
-0.4 
,+0.4 
-0.4 

9/11 



2.7; 

23. 2 1 



0.98 f 

+0.2 
-0.3 
+0.3 
-0.3 
+0.04 
0.04 

43/31 



2.9 
23.1 
1.26 



0.98-f 

+0.1 

-o.i 

+0.1 
-0.1 
+0.08 
0.08 

51/52 



2.9 
23.1 
1.85 



2.3 



0.98-f 

+ 1.0 

.3* 



22.2 



+ 1.0 



1.43 



0.7 
+0.05 



0.04 

21/27 



158/156 



gets. iMateos et al. (2005) found that the distribution of 
QSO spectral slopes can be described as a Gaussian with 
a mean of ax = 0.98 and a standard deviation a a = 0.21. 
From Table [3] we see that when the spectral slopes are 
allowed to vary, the best fit values of ax for RXJ005734 
and RXJ163308 are relatively close to the mean of the 
ax distribution. However, the best fit values of ax for 
RXJ094144, RXJ121803, and RXJ124913 are all unusu- 
ally low, lying within the lowest 3 percent of the distri- 
bution. This is an unlikely circumstance if the intrinsic 
photon indices of our sample were drawn from the same 
distribution as QSOs in general, and suggests that the cold 
photoelectric absorption does not fully account for the un- 
usual X-ray spectral shapes observed in these objects. In 
comparison, the spectral slopes are fixed in the ionized 
absorber model, and so are by design consistent with a 
typical underlying QSO spectrum. 

Moving to aox, we have calculated this quantity for 
each of the model fits and each of the QSOs. In each case 
the restframe 2500 A flux has been estimated from the 
B magnitude listed in lPage. Mittaz fc Carrera" (20011 ). as- 
suming an optical slope of ao = 0.5. The 2 keV flux is 
determined from the continuum model of the X-ray spec- 
tral fit, i.e. it is the intrinsic continuum flux, corrected for 
the modelled X-ray absorption. 

The distribution of aox in nor- 
mal QSOs is known to depend on lumi- 
nosity (|Vignali, Brandt fe Schneider. 2003 ; 

IStrateva et al.. 20051 ) and there is inevitably a selec- 
tion bias depending on whether the sample is X-ray 
or optically selected. The most appropriate reference 
distribution of apx in the literature is that presented in 
Stra teva et al. (2005) for a sample of 155 Sloan Digital 
Sky Survey QSOs with medium-deep ROSAT coverage. 
However, since this distribution is based on an optically- 
selected sample, while our sample is X-ray selected, we 



have constructed a reference distribution of aox for X-ray 
selected QSOs drawn from th e ROSAT Internat ional 
X-ray Optical Survey (RIXOS; iMason et al.. 2000l ). We 
restrict the reference sample to those RIXOS QSOs at 
z > 1, which corresponds to a similar luminosity range 
to our X-ray-absorbed QSO sample, because RIXOS 
has a similar X-ray flux limit to the ROSAT hard 
spectrum survey from which our targets were drawn. We 
calculated aox for the reference sample in an identical 
fashion to our absorbed QSO sample except that we 
assume no intrinsic X-ray absorption and ax ~ 0.98 
for the RIXOS AGN. The reference sample has a mean 
(aox) = 1-49 and a standard deviation aox = 0.16; the 
distribution of aox for the reference sample is shown in 
Fig. [5] For compari son, the SDSS sample presented in 
IStrateva et al. (20051 ) has a very similar distribution of 
aox, with (aox) = 1-48 and aox = 0.18. 

In Table [3] we see that with the exception of 
RXJ124913, the absorbed QSOs have aox values within 
two standard deviations of (aox) for the reference distri- 
bution. For RXJ124913, the level of discrepancy depends 
strongly on the X-ray spectral model. The ionized ab- 
sorber fit gives the most consistent aox value, at 2.3 stan- 
dard deviations from the mean of the reference sample, 
while the fit with a cold absorber and fitted ax gives the 
most discrepant value of aox = 1.96, > 2.9 standard devi- 
ations from the mean of the reference sample. Indeed, this 
value of aox is larger than any ap x in the RIXOS refer - 
ence sample or the SDSS sample of lStrateva et al. (20051 ). 
suggesting that the model with a cold absorber and the 
continuum slope fitted as a free parameter would also 
require a highly unusual spectral energy distribution for 
RXJ124913. While RXJ094144 has rather typical values 
of aox in Table we note that if we correct for the in- 
trinsic extinction equivalent to E(B— V)=0.2, suggested by 
the UV continuum shape (Section l3.1[) . the aox would in- 
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Figure 3. EPIC spectra of the X-ray absorbed QSOs together 
with the ax = 0.98 power law and ionized absorber model 
(black stepped line) and free-a power law and cold absorber 
model (grey stepped line). All spectra are shown in the observer 
frame. Both model and data have been divided by the product 
of the effective area and Galactic transmission as a function of 
energy. 



crease by 0.40, leading to an abnormal aox for the cold 
absorber, fitted ax model for this source as well as for 
RXJ124913. 

Overall, the ionized absorber model is consistent with 
the five QSOs having unremarkable underlying X-ray con- 
tinua and optical to X-ray spectral energy distributions. In 
contrast, if the X-ray absorbers are modeled as cold gas, 
the QSOs require a combination of cold absorption, un- 
usual X-ray continua, and unusual optical to X-ray spec- 
tral energy distributions. We therefore consider that the 
ionized absorber model provides a simpler, less contrived 
solution than the model in which the X-ray absorber is 
cold. We also note that the presence of CIV absorption 
lines with EW> 5 A in the restframe UV spectra of all five 
QSOs (and Si IV in RXJ124913 and RXJ163308) provides 
independent confirmation that these objects are viewed 
through columns of ionized gas. Therefore in what follows 
we adopt the ionized absorber fits as the best description 
of the X-ray spectra. 



4 DISCUSSION 

4.1 Characteristics of the absorbers 

Our study has revealed absorption from ionized gas in 
both the rest-frame UV spectra and in the X-ray spec- 
tra of our sample of X-ray absorbed QSOs. In the 
nearby Univers e, ionized absorbers are common in S eyfert 
galaxies (e.g. iRevnolds. 19971 ; IGeorge et al.. 19981). and 
have been studied in s ome detail l|l31ustin et al.. 20071 ; 
ISteenbrugge et al.. 20051 ). Therefore in trying to under- 
stand the properties of our sample of X-ray absorbed 
QSOs, it is instructive to compare them with those found 
in Seyferts. Amongst Seyfert galaxies, ionized absorption 
is observed primarily in broad-line, (type 1) objects (e.g . 
IRevnolds. 19971 ; IGeorge et al.. 19981 ; IBlustin et al.. 20051 ); 
indicating that the ionized absorbers lie within the ioniza- 
tion cones of the AGN, assuming the standard AGN geo- 
metric unification scheme in which the central regions ar e 
surrounded by a dusty, obscuring torus (|Antonucci. 1 993). 
The presence of ionized absorption within our X-ray ab- 
sorbed QSOs is therefore in keeping with a normal pole-on 
orientation for these objects. 

The ionized column densities, and the degree of 
soft X-ray attenuation, of our X-ray absorbed QSOs 
are considerably higher than are found in the majority 
of nearby Seyferts: with the exception of RXJ163308 
for which the X-ray absorber properties are poorly 
constrained, all of our sample have best fit log Nh 
> 22.5, wh ereas only 3 of the 23 AGN in the sample 
studied by IBlustin et al. (20051 ). have log N H > 22.5. 
Nonetheless, a few nearby objects such as NGC3783 
do have com parable column d en sities and ioniza- 
tion states IBlustin et al.. 20021; iNetzer et al.. 20031: 



iReeves et al.. 20041 ; iKrongold et al.. 20031 ). The ion- 
ized absorbers in Seyferts typically contain ions 
covering a large range of ionization, and where 
good data are available they almost always require 
multiple phases with different ionization param- 
eters to model their spectra (Blustin et al.. 2005 ; 



iMcKernan. Yaqoob fc Reynolds. 20071 ). Moving to higher 
luminosity objects, ionized absorbers appear to be com- 
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Figure 4. Confidence contours for £ and Nh in the ionized absorber model fits to the XMM-Newton spectra. 



mo n in QSOs, with the studies of IPorauet et al (20041 ) 
and iPiconcelli et al (20051 ) suggesting that around 50 per 
cent of low redshift QSOs show some evidence for ionized 
absorption in their X-ray spectra. However, column 
densities of ionized material similar to those found in our 
sample of X-ray abso rbed QSOs are relatively rare, with 
IPiconcelli et al (20051 ) finding that only 2 of their sample 
of 40 nearby QSOs have ionized absorbers with best fit 
log Nh > 22.5. In the small number of nearby QSOs which 
have well-characterised ionized ab sorbers with Nh > 
10 22 cm -2 , notably PG 1114+445 (lAshton et al.. 2004h . 
PDS456 dReeves. O'Brien fc Ward, 20031 ). H0557 - 
385 jAshton et al.. 200rj ). and PG 1211+143 
(|Pounds et al.. 20031 ). the absorbers are found to 
have multiple ionization phases, in common with the 
well-studied Seyfert warm absorbers. 

On this basis, we can expect that the parameters de- 
termined for the single phase absorbers in our sources are 
therefore likely to be representative only of the dominant 
phase of absorption, and the high-£ absorption which we 
have detected is likely to be accompanied by significant 
soft X-ray absorption due to material with log £ < 1.5. In 
fact the multi-ionization-phase nature of the ionized ab- 
sorbers is implied by the absorption lines detected in their 
UV spectra. The best fit ionized X-ray absorber models 
for all 5 objects produce insignificant C IV absorption lines 
(EW < 1 A) , in contrast to the large EW lines (> 5 A) that 
are measured in their rest-frame UV spectra, implying 
that at least two ionization phases of absorber are required 
to produce both UV and X-ray absorption. Note that 
while absorption fr om C IV is seen in about half of QSOs 
(|Vestergaard. 200 3"). CIV absorption lines with equivalent 
widths > 5 A are relat ively rare. Only 3 out of 114 non- 
BAL QSOs studied by I Vestergaard (20031 ) show CIV ab- 
sorption with such large EW. Combined with the statis- 
tic that ~ 15 percent of QSOs have CIV BALs, which 
by definition have CIV EW > 5 A (jGibson et al.. 20091 ; 



in 













n , 
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Figure 5. Reference ctox distribution for z > 1 RIXOS X-ray 
selected QSOs, as described in Section 13.21 

iTolea. Krolik fc Tsvetanov. 2 002) we can conservatively 
estimate the incidence of C IV absorption lines with EW> 
5 A to be 20 per cent. The probability of all 5 of our X- 
ray absorbed QSOs having such absorption lines if they 
were drawn randomly from the population is therefore 
0.2 5 = 3 x 10 -4 , and hence the connection of the CIV 
absorbers with the X-ray absorption is quite secure. 

Our finding that the X-ray absorbers in our sam- 
ple are highly ionized allows us to address the puz- 
zling absence of strong UV attenuation by dust, de- 
spite the large column densities of X-ray absorbing gas. 
Nearby Seyfert 1 galaxies with warm absorbers also show 
this property: the column densities of X-ray absorb- 
ing gas can be large but they have bright ultraviolet 
continua with broad emission lines, e.g. NGC 3783 has 
\o%Nh > 22.5 (jNetzer et al.. 20031 ; [Reeves et al.. 20041 ) as 
well as a strong UV contin uum and broad UV emission 
lines with little reddeni ng (|Crenshaw fc Kraemer. 200ll ; 
lEvans fc Koratkar. 2001 ). 
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As a corollary to the apparent similarity between the 
ionized absorbers found in our X-ray absorbed QSOs and 
some nearby Seyferts, we might wonder whether star for- 
mation is associated with such X-ray absorption in objects 
such as NGC 3783, as it is in X-ray absorbed QSOs. A sys- 
tematic investigation of this issue is beyond the scope of 
this paper, but it is interesting to note that NGC 3783 has 
a 2-10 keV X-ray luminosity of 1.3 x 10 43 ergs s -1 and a n 
infrared luminosity of 3 x 10 10 L (|Mulfanev et al.. 201 J ), 
so that if we scaled its SED to match the X-ray luminosi- 
ties of our X-ray absorbed QSOs its infrared luminosity 
would be in the regime of ultraluminous infrared galaxies, 
commensurate with a substantial star formation rate. 



4.2 Comparison to soft X-ray weak QSOs 

It is interesting to compare the X-ray absorbed QSOs with 
another population of unusual QSOs: those that are classi- 
fied as soft X-ray weak. These objects are characterised by 
large values of a o x , and represent a ro und 10% of optically 
selected QSOs (|Laor et al. 1997m . iBrandt et al. (2000 ) 
found a strong correlation between a ox and the equivalent 
width of C IV absorption, which is found to be > 4A 
in the majori ty of soft X-ray w eak QSOs. This was 
interpreted by IBrandt et al. (2000) as evidence that soft 
X-ray weak QSOs are such because they are absorbed 
in the soft X-ray band. As shown in Section 13.11 the 
QSOs in our sample all have C IV absorption lines with 
EW> 4 A, similar to the CIV lines found in soft X-ray 
weak QSOs. More recently, XMM-Newton observations of 
soft X-ray weak QSOs have indicate d that this absorption 
is primarily due to ionized gas ( Schartel et al.. 2 005; 

Brinkmann. Papadakis fc Ferrero. 20041; 

Piconcelli. Jimenez- Bailon fc Guainazzi. 20041 ). with ion- 
ization parameters ranging from log £ = 1.9 to log £ = 3.2 
and column densities ranging from 3 x 10 22 cm" 2 to 



These column densities and ionization 
parameters are remarkably similar to those of our X-ray 
absorbed QSOs. Given these similarities, we might expect 
that intense star formation would be as common in soft 
X-ray weak QSOs as it is in X-ray absorbed QSOs. To 
our knowledge there has been no systematic study of the 
star formation rates in soft X-ray weak QSOs, but there 
are certainly well known examples of soft X-ray weak 
QSOs which are rapidly forming star s, notably Mrk 231 
and IRAS 07598 + 6508 (jLipari. 19941 ). 



4.3 Location of the absorbers 

We now consider the location, physical nature, and impli- 
cations of the ionized absorbers found in our sample. To 
estimate the location we assume that the absorber is in 
the form of an outflowing wind originating at a distance 
Rabs from the continuum source with density n oc r~ 2 . 
From the definition of the ionization parameter 
£ — L lon /(nr 2 ) (1) 
where Li 0n is the 1-1000 Rydberg ionizing luminosity in 
erg s _1 , such a wind will have an ionization parameter £ 
that is constant with radius r. To take account that the 
winds from AGN are likely to be clumpy or filamentary, we 
characterise the outflow with a macroscopic volume filling 



factor / < 1. The column density of the absorber is dom- 
inated by material close to the base of the outflow, and 
hence R a bs can be taken as the characteristic distance of 
the absorber from the continuum source. Integrating the 
density along the line of sight, we obtain Nb = fRabsn a bs 
where n aos is the density at the base of the outflow. Sub- 
stituting n a bs for n in Equation [T] we obtain the relation 



Rab 



Lion f 



(2) 



To estimate Lion we assume that below 0.35 
keV the SED has the functional form given by 
iMathews fc Ferland ( 19871 ) for a QSO ionizing spec- 
trum, in which the X-ray power law steepens to a = 3 
in the range 56-350 eV, has a slope of a = 1.0 in the 
24-56 eV ra nge, below which it has a s lope of ao = 0.5. 
This SED of IMathews fc Ferland (1987D was constructed 
to be appropriate for objects of QSO luminosity, and 
when joined to our ax =0.98 spectrum at 0.35 keV it has 
aox = 1-47. 

For the fixed ax =0.98 which we have taken in our 
ionized absorber model fits, Li 0n is a factor of 44.3 times 
larger than the restframe 2-10 keV luminosity. For the fill- 
ing factor / we assume a value of log / = —2, which is 
the average log / for ioniz ed absorbers with £ > 0.7 in 
Seyfert galaxies found by lAshton et 80772 006). The de- 
rived values of Li 0n and R a bs are given in Table 2] With 
the exception of RXJ163308, for which R a bs is very poorly 
constrained, the best fit values of R a bs range from 1 to 15 
parsecs. The uncertainties given for R a bs are derived from 
the statistical uncertainties on the product £ x Nh, and so 
do not include the uncertainty on the volume filling fac- 
tor. As an estimate of thi s additional uncertainty, we note 
that lAshton et al. (2006) obtained a standard deviation 
in log / of 0.8 for Seyfert warm absorbers with £ > 0.7. 
For comparison, we also provide in Table [4] the locations 
of the inner edge of the dusty toru s Rtorus for our s ources, 
calculated using equation 5 from iBarvainis" (1987), with 
a dust sublimation temperature of 1260 K (as implied 
by the composite near-IR QSO spectrum obtained by 
iGlikman. Helfand fc White. 20061 ) and taking the ultra- 
violet luminosity to be equal to Li 0n . We see that R a bs is 
in all cases smaller than, or consistent with, Rtorus- 

It would therefore be natural to assume that the ab- 
sorbers originate in the AGN themselves, rather than in 
their host galaxies. This solution is attractive, because it 
is compatible with the lack of optical extinction in these 
objects: if the absorber is driven as a wind, either from 
the accretion disc or from evaporation of the inner edge 
of the molecular torus, then dust will be sublimated be- 
fore (or as) it enters the flow. With this in mind, it is 
instructive to obtain upper limits for the distance of the 
absorbing gas from the nucleus at which it could remain 
ionized by the AGN radiation field, using Equation [2] and 
taking the limit case of filling factor / = 1. Doing so, we 
find that the X-ray absorbers in RXJ094144, RXJ121803 
and RXJ124913 must lie within a few hundred parsecs of 
the ionizing source, within or close to the AGN rather than 
being distributed throughout the host galaxy. 

A key indicator of the importance of an ionized AGN 
wind is the ratio of the mass outflow rate M ov .t to the 
accretion rate M acc . To determine this ratio, we first take 
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Table 4. Characteristic scales and outflow rates for the sample 
of X-ray absorbed QSOs. Li on is the 1-1000 Rydberg ionizing 
luminosity in erg s -1 , R a bs is the distance of the base of the 
absorber from the continuum source in cm, Rtorus is the dis- 
tance of the inner edge of the dusty torus from the continuum 
source in cm, M acc is the mass accretion rate and M ou t is the 
mass outflow rate in the ionized absorber. An '*' indicates a 
confidence limit that is not constrained by the data. Note that 
the uncertainties are purely statistical, and do not include the 
systematic uncertainty on the filling factor /, on which both 
Rabs an( l M ou t depend linearly. As explained in Section 1 1.31 
we estimate an uncertainty on log / of ~ 0.8. 



Object 

RXJ005734 
RXJ094144 
RXJ121803 
RXJ124913 
RXJ163308 



logL ion logR abs logR torua log(M out /M acc ) 



46.7 
46.2 
46.6 
46.8 
47.0 



19 6 +0 - 7 

18 4 +0 - 6 
+0.5 
0.6 
+0.3 
0.2 

2.4 



18.6 

18.8_ 

20.7^ 



19.4 
19.2 
19.4 
19.5 
19.5 



1.1 



+0.4 
0.3 



-0 1+ 04 

9+ 02 

1 2+* 



Lion to be 0.6 x L acc (jVasudevan fc Fabian. 200*91 ) where 
Lace is the bolometric accretion luminosity, so that 

Lion — 0.6Lacc — 0.6eMaccC 2 (3) 

where c is the speed of light, and e is the accretion effi- 
ciency. We can obtain the mass outflow rat e of the wind 
using equation 18 from lBlustin et al. (20051 ): 

1 .23mp L i on f vQ 



Man 



(4) 



where m p is the rest mass of a proton and v is the outflow 
velocity, and is the solid angle of the outflow. Combining 
equations [3] and [4] we obtain: 

Mont/Mace = 0.74m p fvnec 2 /t (5) 
We take e = 0.1 and as before we assume / = 0.01. 
We take the opening angle of the outflow Q = n stera- 
dians, corresponding to the opening angle of the torus in- 
ferred fr om the type 1 / type 2 number ratio in nearby 
Seyferts ([Maiolino fc Rieke. 1 995). For the velocity v we 
assume that the X-ray absorbing ionized outflow is asso- 
ciated with the fastest outflow identified in the restframe 
ultraviolet (Table [2}. This is a reasonable estimate, since 
the highest ionization absorbers are typically (though not 
always) the hi ghest velocity sys t ems in individual Seyfert 
galaxies (e.g. Sako et al.. 20031 : ISteenbrugge et aL , 2005; 
ICostantini et al., 2007^ At the same time it is a conser- 
vative estimate (i.e. may underestimate the mass outflow 
rate) for two reasons: firstly, the signal to noise ratios 
of our optical spectra (particularly for RXJ094144 and 
RXJ121803) are not high enough that we can be confi- 
dent we have identified all of the absorption line systems, 
and secondly our estimates of the systemic velocities are 
based on the UV broad emission lines, which themselves 
are frequently obser ved to have a net outflow velocity of 
order 1000 km s" 1 (|Espev et al.. 19891 ). The derived val- 
ues of log Mout/ Macc are listed in Tabled In the case of 
RXJ094144, the C IV absorption line that is identified in 
the rest-frame UV spectrum is redshifted with respect to 
the UV emission lines, so its true velocity with respect to 
the system is too uncertain for us to attempt to estimate 
the mass outflow rate. The mass outflow rates derived 
for RXJ005734, RXJ124913 and RXJ163308 are about 10 
times their mass accretion rates, while in RXJ121803 the 
accretion and outflow rates are comparable. 



4.4 Implications for QSO evolution 

At this point it is worth considering how these re- 
sults bear on the broader picture of how X-ray ab- 
sorbed QSOs relate to the evolu tion of massive b lack 
holes and their host galaxies. In IPage et al~ (20041) we 
argued that the submillimetre properties of X-ray ab- 
sorbed QSOs imply an evolutionary sequence in which 
submillimetre galaxies evolve into QSOs, with the X-ray 
absorbed QSOs representing a brief, transitional phase 
between these two stages. The discovery that submil- 
limetre galaxies typically contain more heavily absorbed, 
lower luminosity AGN than X-r ay absorbed QSOs, sup - 
po rts this evolutionary picture (I Alexander et al.. 20051 ). 
In IStevens et al. (20051 ) we further proposed that the X- 
ray absorbed QSOs are in the process of driving the in- 
terstellar media out of their host galaxies via radiatively 
driven winds, in a ccord with th e models developed by 
ISilk fc Rees (1998TI. iFabian (19991), iGranato et al. f2004f) 
and Di Matteo. Springel fc Hernquist (20051 ). Now, with 
our XMM-Newton observations we have shown that the 
X-ray absorption in these objects is indeed due to sub- 
stantial, ionized, QSO-driven winds. The space density 
of X-ray absorb ed QSOs is ~ 15 per cent that of un- 
absorbed QSOs jPage et al.. 200j ; ISilverman et al.. 20051 : 
IPage et al.. 20061 ) . implying that if all QSOs go through 
an X-ray absorbed phase in their evolution, this phase 
lasts only ~ 15 per cent as long as the unabsorbed phase. 
These relative lifetimes may be a natural consequence of 
the high mass outflow rates in the X-ray absorbed QSOs: 
for the majority of our X-ray absorbed QSOs we estimate 
a mass outflow rate which is ~ 10 times higher than the 
accretion rate, implying that the fuel supply will be de- 
pleted 10 times faster in an X-ray absorbed QSO than in 
an unabsorbed QSO for the same accretion rate and fuel 
reservoir. However, this argument is more complicated if 
the two types of QSO fit within the evolutionary sequence 
discussed above, because in this case the fuel supply for an 
unabsorbed QSO will consist only of the fuel that remains 
after the X-ray absorbed phase. Nonetheless, if material is 
ejected at ~ 10 times the accretion rate for ~ 15 per cent 
of the (X-ray absorbed + unabsorbed) QSO lifetime, the 
total mass of material ejected by the QSO will be of the 
same order as that accreted. 

The origin of the ionized, absorbing gas is then a fun- 
damental issue for the physics of X-ray absorbed QSOs 
and their evolution. Assuming that X-ray absorbed QSOs 
have the usual AGN structure envisaged in geometric uni- 
fication schemes, the two obvious sources of material are 
the accretion disc and the dusty torus. If the wind is 
driven from the accretion disc, the ejected material must 
pass through the accretion disc during the X-ray absorbed 
phase; in this case it is difficult to understand why this 
happens on a timescale which is short compared to the 
lifetime of the QSO. On the other hand, if the ionized 
wind is driven from the torus then the timescale for the 
X-ray absorbed phase corresponds to the time required to 
erode the inner part of the torus a nd so enlarge substan - 
tially its opening angle; as shown bv lKrolik fc Kriss (2 001) 
this timescale can be much shorter than the lifetime of the 
QSO. For accretion at the Eddington rate, the escape ve- 
locities at the distances of the inner edge of the torus (as 
given in Table [4| are between 900 and 1500 kms -1 , and 
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scale inversely with the square root of the accretion rate in 
Eddington units. Hence, of our X-ray absorbed QSOs, only 
the BALQSO RXJ124913 has outflow velocities which are 
too high to plausibly be associated with an outflow from 
the inner edge of the torus, implying an accretion disc 
origin for its outflow. 

Viable AGN feedback models which reproduce 
the M — a relation generally incorporate the feed- 
back in the form of kinetic energy given to gas 
within the imme diate (~100 pc) environment of 
th e AGN (e.g. [Fabian. 19991; iGranato et al., 20041 : 



Pi Matteo. Springe! fe Hcrnquist, 2005; 
Hopkins et al.. 20061 ; ISiiacki et al.. 20071) . which then 
interacts with gas on a larger scale to terminate star 
formation in the host spheroid. It is generally agreed that 
of order 5 per cent of the AGN radiative output must 
be fed back as kinetic energy in order for the models 



to succeed (ISilk fe Rees. 19981; IWvithe fe Loeb. 20031 : 



|Pi Matteo. Springel fc Hernquist. 20051 ). We can make 
a simple estimate for the kinetic energy of the outflows 
compared to the radiative output of the X-ray absorbed 
QSOs as follows. We take the total energy radiated during 
the X-ray absorbed phase to be Et = tM acc c? where e is 
the radiative efficiency, c is the speed of light and M acc is 
the mass accreted during the X-ray absorbed phase. The 
kinetic energy built up by the outflow is Eo = 0.5M out v 2 
where M out is the mass supplied to the outflow during 
the X-ray absorbed phase and v is the velocity of the 
outflow. The fraction of the radiated output of the QSO 
which is fed back to the surrounding gas as kinetic energy 
of the outflow is then 
E _ 0.5 Mout v 2 

~E~r ~ ~~r x Ma^r x ~& [b) 

If we take M out /M acc =M out /M acc ~ 10 for the X-ray 
absorbed phase (as we find for the majority of our X-ray 
absorbed QSOs), and a typical outflow velocity of 8000 
km s _1 (see Table[2}, we obtain Eo/ Et ~ 4 per cent. The 
value of Eo/ Et so derived is independent of the radiative 
efficiency e, because the e in Equation HJ] cancels with that 
in Equation [5] from which M ou t/M aC c is obtained. On the 
other hand, the uncertainty on M ou t/ M aC c is in truth very 
large (> a factor of 10), depending on the unknown filling 
factor /, and the square of the outflow velocity v, which 
may not be the same for the X-ray absorbing phase as the 
UV absorber. Nonetheless, Eo/Et ~ 4 per cent is quite 
consistent with the level of feedback required to terminate 
star formation and produce the M — a relation, implying 
that the ionized outflows in X-ray absorbed QSOs could 
plausibly be the mechanism through which star formation 
is terminated in massive galaxies. 



5 CONCLUSIONS 

We have presented X-ray spectra from XMM-Newton 
EPIC and rest-frame ultraviolet spectra from ground- 
based telescopes for five X-ray absorbed, submillimetre- 
luminous QSOs. All five QSOs exhibit strong C IV ab- 
sorption lines in their ultraviolet spectra with equivalent 
width > 5 A. The X -ray spectra can be modelled success- 
fully in terms of x 2 with either cold or ionized absorbers. 
The cold X-ray absorber model requires that the QSOs 



have unusually flat X-ray continuum shapes and unusual 
optical to X-ray spectral energy distributions, while the 
ionized absorber model does not require abnormal under- 
lying continuum properties. This finding, coupled with the 
presence of strong C IV absorption lines in the UV leads 
us to favour the ionized absorber model over the cold ab- 
sorber model. Assuming that the X-ray absorbing gas is 
outflowing with the same velocities as the C IV absorbers, 
we are able to investigate the likely location, mass outflow 
rates and energetics of the ionized absorbers. We find that 
the X-ray absorbing gas is likely to be located within 10 pc 
of the continuum source, and so is associated with the 
active nucleus rather than the surrounding host galaxy's 
interstellar medium. We estimate that the fraction of radi- 
ated power that is converted into kinetic luminosity of the 
outflowing winds is typically ~ 4 per cent, in agreement 
with estimates for the kinetic feedback from QSOs that 
is required to produce the M — a relation. This finding is 
thus consistent with the hypothesis that X-ray absorbed 
QSOs represent the transition phase between obscured ac- 
cretion and the luminous QSO phase in the evolution of 
massive galaxies. 
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